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Summary 

A procedure is described for designing dual- 
reflector antennas. The analysis is developed by tak- 
ing each reflector to be the envelope of its tangent 
planes. The slopes of the emitted rays are speci- 
fied rather than the phase distribution in the emit- 
ted beam. Thus, both the output wave shape and 
the angular distribution of intensity can be specified. 

Computed examples include variations from both 
Cassegrain and Gregorian systems. These examples 
include deviation from uniform source distributions 
and from the parallel-beam property of conventional 
systems. 

Introduction 

In theory it is possible to specify, within lim- 
its, both the emitted amplitude distribution and the 
phase distribution of a dual-reflector system when 
both surfaces are properly shaped. The problem of 
determining these shapes has been treated in refer- 
ences 1 and 2 with a combination of differential and 
algebraic equations. This system of equations tends 
to be unwieldy, and consequently in both references 
the method of solution is only indicated, with the 
specific formulas omitted. 

Although the present approach is an approximate 
procedure, it can, in theory, be made as accurate 
as desired by taking sufficiently small step sizes. 
The simplifying concept is to treat each of the two 
reflectors as the envelope of its tangent planes. This 
technique permits the mathematical problem to be 
reduced to one of solving a set of nonlinear algebraic 
equations. 

Computed examples include both modified 
Cassegrain and modified Gregorian systems. Inas- 
much as the emphasis in the present analysis is on 
simplicity of concept, only axisymmetric systems are 
treated. It should be noted that if the reflector sys- 
tem utilizes only a segment (e.g., a quadrant) of the 
axisymmetric design, then a cross polarization exists 
in the output beam since the antenna normally does 
not emit a circularly polarized wave. 

Symbols 

b exponent of cos 6q 

E beam energy 

/ intensity 

m slope of ray relative to 

system centerline 

_ -2m 2 
— m%-l 


x,y 

coordinates, x taken along 
system centerline and y in 
radial direction 

(X l ,Y l ),(X 2 ,Y 2 ) 

midpoint of straight seg- 
ment tangent to reflector 
meridian line 

x s ,Y 3 

point at which ray intersects 
reference plane 

{xi,yi),(x 2 ,y 2 ) 

initial point of straight 
segment tangent to reflector 
meridian line 

e 

angle ray makes with 
system centerline 

Subscripts: 

0 

origin or ray emanating 
from origin 

1 

tangent segment of 
subreflector 

2 

tangent segment of main 
reflector 

3 

vertical plane along which 
output intensity is specified 

b 

ray reflected from main 
reflector 

c 

ray reflected from subreflec- 
tor to main reflector 

i 

index 

min, max, total 

minimum, maximum, and 
total 


Analysis 

Reflector and Ray Geometry 

As was mentioned in the Introduction , each 
curved reflector surface is treated as the envelope 
of its tangent planes. Since each of these surfaces 
is axisymmetric, it can be specified by the meridian 
line cut by a vertical plane. Then the tangent plane 
along this meridian appears simply as a straight line 
segment. (See fig. 1(a).) 

The geometry for tracing a ray through the sys- 
tem is shown in figure 1(b). The reflection condition 
at the subreflector gives 


R 


2 6\ — 0 C + (7r + $o) 


(i) 



(b) Angles and slopes. 
Figure 1. Basic system geometry. 



To find the relation between the slopes, we take the 
tangent of both sides: 


tan 20i = 


2 tan 6\ 


tan 0 C + tan ( 


1 — tan 2 1 — tan 0 C tan 0q 
Similarly, at the second reflector, 

tan( 2 <? 2 ) = tan[03 + ( 0 C - rr)] 


or 


2 tan 02 tan 0 3 + tan 0 C 


( 2 ) 


(3) 


(4) 


1 — tan 2 02 1 — tan 03 tan 0 C 

Since the slope of the local tangent to the reflector 
is the negative reciprocal of the slope of its normal 
(mi = — cot#i), equation (2) yields 


-2mi mo + m c 


to 2 — 1 


1 — mom c 

and, similarly, equation (4) yields 

-2m2 _ m3 + m c 
1 - m3m c 


m\ - 1 


(5) 


( 6 ) 


Referring to figure 1(a), if we denote the coordi- 
nates of the initial point of a subreflector segment by 
(ri, yi) and denote the reflection points by capital 
letters (Xi, Yi), then the reflection point is specified 
to be at the midpoint of the segment simply by tak- 
ing the segment length to be twice the distance from 
(zi)2/l) to (Xi,Yi). The main reflector is treated 
similarly. The slopes and point coordinates are re- 
lated linearly as follows: 


m Y2 — Y\ 

c X 2 - Xi 

(7) 

™ Yi~yo 

m ° V 

X\ - Zo 

( 8 ) 

™ - Y1 ~ yl 
m l - V 

Xi - Zi 

(9) 

_ Y2 - V 2 

m 2 V 

A 2 - X2 

(10) 


The intensity distributions are determined as fol- 
lows. The intensity distribution to be assigned ar- 
bitrarily as the output of the system is most conve- 
niently specified along a vertical plane, usually taken 
near the aperture plane. Thus, if the output inten- 
sity distribution /(Y3) is assigned as a function of 
Y3 along a vertical plane taken at X3, then the ray 
that intersects this plane at Y3 is related to the ini- 
tial source emission angle 9 q through the energy rela- 
tion, which is given in the next section. Furthermore, 


it is appropriate, with the present method, to spec- 
ify the ray slope distribution of the system output 
beam rather than the phase distribution. Thus (see 
fig. 2(a)), in the equation for the ray reflected from 
(X 2 ,Y 2 ) through (X 3 ,Y 3 ), 


m b = 


Y3-Y2 


X3-X2 

and X3, Y3, and m b are all known quantities. 


(11) 


Energy Relation 

The value of I3 is determined from the energy 
relation as follows. The intensity distribution Io(Oq) 
emitted by the source antenna is (see fig. 2(b)) 

2ir I o( 6 q) sin Oq (16q 

and consequently all the energy emitted within this 
segment is 

fO 0 

£o(0o) = 2 tt / I o {0) sm6 d6 (12) 


and the relative amount of energy emitted is 

Wo) _ Wo) 

Eq (00 ,max) -^total 


(13) 


where 0o,max denotes the edge of that part of the 
source beam that is to be utilized. 

The intensity of the output beam is specified as a 
function of Y3 along the vertical plane at X3. Thus, 
the energy emitted through the annulus at Y3 with 
width (IY3 is (see fig. 2(a)) 


27r/3(Y 3 )Y 3 dY z 

The energy emitted through the ring at Y3 is 

fY 3 

E(Y 3 ) = 2tt / 

JY'. 


’Y 3 ,z 


h(y)y dy 


(14) 


where the ray through Y3 )in j n is that emitted at 
#0,min by the source. The relative energy is 


E(Y 3 ) _ E(Y Z ) 


E(Y Z ,max) -^total 


(15) 


Comparing equation (15) with equation (13) enables 
one to determine the Y3 corresponding to a given Oq. 
As a rule, the prescribed intensity distribution I Z (y) 
in equation (14) is taken as a relatively simple an- 
alytic expression. Consequently, the energy integral 
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(a) Vertical reference plane at X3. 



(b) Source rays at unit circle. 

Ray (constructed 
normal to wave front) — v . 


Prescribed 
wave shape 
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(c) Method for determining Y3 and Of, from prescribed wave shape. 

Figure 2. Geometry for identifying source ray at Oq with emitted ray through (X3, I 3 ). 


in equation (14) can be evaluated in closed form and 
E(Y 3 ) is determined in analytic form. Furthermore, 
the source intensity distribution I${6) can often be 
described as some power of cos 6 or as a linear com- 
bination of such functions, so that E(0$) can be de- 
termined in analytic form from equation (12). 

Before proceeding to the solution of the equations 
for the reflector surfaces, we may pursue further the 
significance of the form of the system performance 
specification. Although for some problems it is 
appropriate to specify the intensity and beam 
direction distributions along some vertical plane, it is 
important to observe that these distributions can be 
obtained by specifying the more fundamental quan- 
tities of output wave shape and the intensity dis- 
tribution as a function of ray direction. Thus, in 
reference to figure 2(c), the wave shape and the rela- 
tive energy as a function of 6 b are prescribed. Since 
the wave shape is known, normals to the wave shape 
can be constructed and their intersections with the 
vertical plane at X 3 determined. These normals, 
which represent rays, have slope m b = tan0&. Thus, 
with E(0 b ) prescribed, these intersections determine 
E(Y 3 ). Consequently, both m b and Y3 are deter- 
mined for each ray of the output beam. 


Solution of Equations 

To determine the reflector shapes, the set of seven 
equations (5) to (11) are to be solved for the unknown 
quantities X\,Y \, and m c . The ap- 
proach taken here is to eliminate all unknowns except 
m2, solve the resulting (highly nonlinear) equation 
for m 2 numerically, and then substitute back into 
the other equations to determine the remaining un- 
knowns. The details of this procedure follow. 

Define 

R{m 2 ) = 2 2?W2 (16) 

77^2 1 

and substitute R into equation (6), which (noting 
that m3 = m b ) can then be written as 

(1 - m b m c )R{m 2 ) = m b + m c 


and solved for m c to obtain 


m c (m 2 ) = 


R(m 2) - m b 
m b R(m 2 ) + 1 


(17) 


Equations (10) and (11) are each solved for Y 2 , 
and then Y 2 is eliminated by equating the resulting 
expressions: 


This equation is solved for X 2 to yield 

= - n - + u (18) 

m 2 - m b 

This expression may be substituted back into equa- 
tion (11) to obtain an expression for >2 as a function 
of m2 only: 

*2(^2) = Y 3 - m b [X 3 - X 2 (m 2 )] (19) 

Eliminating Y\ between equations (7) and (8) yields 
Y 2 - m c (X 2 - X\) = yo + m 0 (Xi - x 0 ) 


which can be solved for X\ and expressed as a func- 
tion of m2 with substitutions from equations (17), 
(18), and (19): 

V X m c (m 2 )X 2 (m 2 ) - Y 2 (m 2 ) - m 0 x 0 + y 0 

Xlim2) = m„(m 2 )-m 0 

( 20 ) 

This result is substituted back into equation (8) to 
obtain 


Yi(m 2) = Vo + m 0 [Xi (m 2 ) - zo] 
Equation (9) now becomes 


mi(m 2 ) — 


U(m 2 ) - yi 
Xi (m 2 ) - x\ 


( 21 ) 


( 22 ) 


Substituting from equations (17) and (22) into equa- 
tion (5) yields 


- 2mi(m 2 )[l - mom c (m2)] + [mo + m c (m 2 )] 
x {l - [mi(m 2 )] 2 } = 0 (23) 

This equation is solved numerically by a forward 
seeker algorithm that finds the zeros of the function 
on the left side of equation (23). Once the root is 
found, we can find the remaining quantities by re- 
peating the above substitutions with the known value 
of m 2 . Thus, m c is obtained from equation (17), and 
X 2 ,Y 2 , X i, and Y\ are obtained from equations (18) 
to (21). 

To determine the next pair of points on the two 
reflectors, Oq is incremented. Then xi, y\, x 2 , and 
y 2 are incremented by specifying (X\, Y,) to be the 
midpoint of a segment on the subreflector and sim- 
ilarly for ( X 2 ,Y 2 ) on the main reflector. Thus, for 
example, 


*3 ~ (^3 ~ X 2 )m b = y 2 + m 2 (X 2 - x 2 ) 


*1 ,* = + 2 (^l, i-1 ~ 
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The procedure can then be repeated for the new 
value of Oq. Inasmuch as equation (23) is highly 
nonlinear and possesses multiple roots, some care 
must be exercised in setting the limits of the range 
of m 2 over which the numerical algorithm seeks a 
solution. Fortunately, the limits are known to a 
close approximation because specifying even sizeable 
variations from a uniform intensity distribution or 
from a parallel-ray output beam does not result 
in large geometry variations from a conventional 
Cassegrain or Gregorian system. 

Computed Examples 

Figure 3 shows a modified Cassegrain system for 
which the source distribution varies as cos 8 6q but 
the emitted beam is specified to have a uniform 
distribution. Figure 4 gives a similarly designed 
modification of a Gregorian system. 

Figure 5 shows a modified Cassegrain system for 
which the slopes of the rays of the output beam are 
specified to increase gradually up to a value of 0.25: 

m 3 = 0.25 (V^-) 

\ ^3, max J 

Figure 6 shows an offset system consisting of seg- 
ments of a subreflector and a main reflector designed 
so that the emitted rays all have a slope of 0.25 rel- 
ative to the system centerline. In such a system, the 


cross-polarization phenomenon mentioned in the In- 
troduction would exist. 

Concluding Remarks 

A procedure has been described for designing 
dual-reflector antennas. The analysis was developed 
by taking each reflector to be the envelope of its 
tangent planes, so that the reflection condition is 
satisfied on each of these planes. The slopes of the 
emitted rays were specified rather than the phase 
distribution in the emitted beam. 

Computed examples included variations from 
both Cassegrain and Gregorian systems. These ex- 
amples include deviation from uniform source distri- 
bution and from the parallel-beam property of con- 
ventional systems. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 2, 1988 
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Figure 5. Modified Cassegrain system with ray slopes of output beam specified to increase to 0.25. 
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